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基于 改进 匈牙利 算法 的 航路 网 络 延误 优化 


王 莉莉 ， 王 航 臣 
(中 国民 航 大 学 天 津 市 空 管 运行 规划 与 安全 技术 重点 实验 室 ， 天津 300300) 


摘 要 : 针对 我 国航 路 网 络 流量 分 配 依靠 管制 员 经 验 ， 缺 乏 量 化 模型 协助 决策 这 一 问题 ， 对 航路 网 络 的 延误 优化 模 
型 及 其 求解 算法 展开 研究 。 首 先 ， 建 立 了 目标 函数 为 总 航行 时 间 最 小 、 带 容量 约束 的 流量 指派 模型 ; 其 次 ， 针 对 经 
典 匈牙利 算法 只 适用 于 中 小 型 运算 的 问题 ， 在 效率 矩阵 的 计算 上 作出 加 图 的 零 元 素 个 数 最 大 概率 与 效率 和 矩阵 的 阶 数 
相等 的 改进 ; 再 次 ， 根 据 航路 运行 效率 随时 间 和 流量 变化 的 特点 ， 使 用 MVMVC 排队 论 模 型 构建 了 航路 的 费用 函数 ， 
用 变化 的 费用 替代 效率 矩阵 中 不 变 的 效率 值 ， 使 效率 和 天 阵 随 时 间 和 流量 变化 ; 最 后 ， 以 中 南 空 管 局 与 西南 空 管 局 的 
部 分 空域 2016 年 11 月 的 实际 运行 数据 作为 算 例 。 结 果 表 明 ， 在 通行 能 力 优化 方面 ， 算 法 能 将 进 场 通行 能 力 提升 
8.372%， 离 场 通行 能 力 提升 8.999%; 在 延误 优化 方面 ， 能 降低 每 架 航 空 器 平均 延误 ; 在 算法 性 能 方面 ， 具 有 比 经 典 
匈牙利 算法 迁 代 次 数 少 、 求 解 时 间 短 的 优点 ， 更 适用 于 实际 的 管制 运行 。 
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Optimization of en-route network delay based on improved Hungarian algorithm 


Wang Lili, Wang Hangchen 
(China Tianjin Key Laboratory of ATC Operation Planning & Safety Technology, Civil Aviation University of China, Tianjin 
300300, China) 


Abstract: Due to cope with relying on controllers badly and lacking of quantitative models to assist decision-making in air 


= traffic flow management, the delay optimization model and its solving algorithm for route network are studied. Firstly, a 
traffic flow assignment model with the objective function of minimum flight time and capacity constraints was established. 
Secondly, classical Hungarian algorithm is only suitable for small and medium-sized operations, the improvement of the 
efficiency matrix was made to improve the algorithm that the maximum probability of the number of zero elements was 
equal to the order of the efficiency matrix; Next, the cost function of the en-route was defined by combining the M/M/C 
queuing theory model. The function, which replaces the constant efficiency value in the efficiency matrix with a variable 
cost, made the efficiency matrix change with time and flow. Finally, part airspaces of China and the actual operational data 
of 2016 were taken as examples. The results show that the algorithm can increase the access capacity by 8.372% and the 
departure capacity by 8.999%; In terms of delay optimization, the average delay of each aircraft can be reduced; in terms of 
algorithm performance, it has less iterations in advantages and less solution time than the classical algorithm, and is more 
suitable for actual control operation. 


Key words: air transportation; delay optimization; Hungarian algorithm; en-route network; assignment problem; cost 


function 

0 引言 行 建 模 ， 但 只 能 精确 地 求解 小 规模 的 航路 网 络 品 , 对 于 大 规模 
六 的 问题 ， 只 能 采用 启发 式 算 法 B41 或 将 问题 分 解 简化 577， 从 
航路 网 络 作 为 空中 交通 系统 的 运行 载体 ， 其 运行 效率 的 而 减少 问题 的 运算 时 间 。 航 路 网 络 的 运行 优化 ， 无 论 是 以 地 
高 低 直 接 决 定 了 系统 的 运行 效率 山 。 管 制 员 依靠 自身 经 验 管 等 待 时间 最 短 !9、 空 中 等 待 时 间 最 短 !69 还 是 航行 时 间 最 
理 空中 流量 ， 缺 乏 量化 的 模型 为 管制 员 向 各 航路 分 配 交 通 流 短 00、 飞 行 油耗 与 排放 最 低 404， 本 质 都 是 有 限 的 时 间 与 空 牛 
提供 决策 依据 ， 这 在 大 流量 的 复杂 空域 易 造 成 管制 员 疲 劳 ， 资源 指派 与 分 配 的 组 合 优化 问题 。 虽 然 如 今 计 算 机 的 运算 能 
听 待 一 个 能 描述 空域 特点 的 方法 协助 管制 员 决 策 。 空 中 交通 力 已 有 大 幅 提 升 , 但 大 规模 ATFM 的 组 合 规模 使 得 在 实际 管 

流量 管理 (air traffic flow management, ATFM) 问题 是 空中 交 制 运行 中 ， 很 多 算法 无 法 投入 应 用 。 
通 管制 中 的 经 典 问题 ， 对 此 各 国 专 家 学 者 建立 了 各 种 目标 函 本 文 针 对 网 络 流 算 法 解决 ATFM 时 间 复 杂 度 高 {9 与 经 典 
数 和 约束 条 件 不 同 的 优化 模型 ， 并 根据 所 建 模型 的 特点 ， 设 匈牙利 算法 和 欠 代 次 数 多 ， 只 适用 于 中 小 型 指派 问题 5 的 缺 
计 不 同 的 算法 求解 。 点 ，a) 建 立 了 一 个 带 容 量 约束 的 0-1 整数 规划 模型 ， 并 将 其 
现 有 的 短期 ATFM 的 研究 ， 多 采用 数学 规划 方法 或 计算 转换 为 指派 问题 求解 ，b) 在 匈牙利 法 需要 使 用 的 效率 矩阵 的 
机 仿真 的 方法 。 其 中 数学 规划 方法 多 以 网 络 流 理论 为 基础 进 ”计算 上 ， 作 出 了 加 圈 的 零 元 素 个 数 最 大 概率 与 效率 矩阵 的 阶 
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数 相 等 的 改进 ， 从 而 减 小 效率 和 矩 阵 的 迭代 次 数 ，c) 结 合 航 路 
网 络 航路 运行 效率 随 流 量变 化 的 特点 ， 提 出 使 用 随 流量 等 因 
素 变 化 的 阻抗 函数 奉 代 经 典 方法 中 效率 矩阵 中 不 变 的 系数 。 

最 后 通过 中 南 和 西南 空 管 局 所 辖 的 部 分 空域 与 2016 年 11 月 
的 实际 数据 作为 算 例 验证 了 模型 与 算法 的 效率 与 运算 质量 。 
1 ”问题 描述 

1.1 航路 网 络 的 数学 描述 


点 t=0,4.…,k 的 航路 网 络 中 ， 其 中 大 和 大 为 节点 数量 ，|1 
于 我 国 高 度 层 分 配 的 原则 是 “东单 西 双 ” 可 将 其 视 为 平面 有 
向 赋 权 连通 图 。 设 有 向 图 G=(V 包 EC(a)), 记 V={V,Vy} ， 其 


中 Vi 表示 航路 交叉 口 的 集合 ，V 表示 起 旋 点 机 场 的 集合 ; 


上 ={e,,,} 表示 航路 路 段 的 集合 ， 其 中 m 表示 航路 的 起 点 ，n 


表示 航路 的 终点 。 并 使 用 权 和 矩阵 法 描述 航路 网 络 ， 从 而 方便 
后 文 指派 问题 的 计算 机 求解 ， 转 换 方法 可 参见 文献 [13]。 
1.2 通行 能 力 及 其 影响 因素 
航路 网 络 的 基本 通行 能 力 这 里 界定 为 在 既定 的 网 络 结构 
和 管制 运行 规则 下 ， 一 定 交 通 流 到 达 情 况 与 分 布 特征 下 ， 管 
员工 作 负 和 荷 可 承担 的 最 大 通行 量 。 通 行 能 力 的 影响 因素 很 
多 ， 如 天 气 因素 、 军 航 活 动 、 地 形 因素 、 交 叉 点 航线 冲突 、 
限制 区 及 其 他 边界 条 件 限 制 等 03。 引 入 一 个 决策 变量 % ， 其 
取 值 只 能 是 1 或 0。 并 令 
5 i { 通 过 航路 rs 0) 
0, 当 不 指派 航班 通过 航路 rs 
以 往 关 于 ATFM 的 研究 中 ， 常 将 航行 费用 视 为 常数 
168,101， 这 对 于 非 拥挤 的 网 络 流量 管理 是 比较 合适 的 ， 但 对 于 
内 部 拥挤 的 网 络 ， 该 方法 计算 的 结果 与 实际 情况 存在 较 大 出 
入 。 实 际 的 航路 网 络 存在 着 严重 的 繁忙 拥挤 现象 ， 费 用 随 着 
流量 的 增加 而 递增 ， 故 需 设计 一 个 符合 航路 网 络 运 行 特点 的 
一 般 认为 ， 众 多 因素 中 与 航空 器 分 配 关系 最 密切 的 是 时 
间 (是 否 能 降低 航班 的 延误 )、 距 离 (是否 能 提升 航路 网 络 的 
运行 效率 ) 与 相关 费用 (是 否 能 达到 效益 最 大 化 )。 不 同 的 飞 
行 类 型 对 影响 因素 的 关注 程度 有 所 不 同 ， 故 定义 费用 函数 为 
c(fs(D) = out fs (DD) + W,)S (2) 
其 中 : a 表示 流量 系数 ， 代 表 流 量 带 来 的 阻抗 影响 ， 可 由 实际 
数据 标定 ;4 表示 节点 的 服务 率 ;， 大 GD) 表示 终端 区 上 时 刻 从 
节点 7 到 节点 s 的 空中 交通 流量 ; 8 为 调整 系数 ， 可 根据 实 
际 数据 进行 标定 ; xA 表示 航空 器 飞 过 航线 rs 所 需 的 时 间 ， 设 
心 表示 第 rs 段 航线 的 长 度 ，， 表 示 空 中 交通 流 的 平均 速度 ， 
可 表示 为 


法 


十 
i 


E: 
外 


7 


(G3) 


在 确定 了 决策 变量 克 和 各 航路 费用 函数 < 这些 基础 条 件 
后 ， 可 以 开始 考虑 整个 航路 网 络 的 通行 能 力 模型 。 有 文献 证 
明 04 单 位 时 间 内 航空 器 的 到 达 服 从 泊 松 分 布 ， 机 场 节点 的 服 
务 服从 负 指 数 分 布 ， 故 平均 延误 屿 可 用 M/M/e 排队 论 模 型 
中 的 Little 公式 05 求 得 ; 

Ww, = (4) 


= I 
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2 航路 网 络 通行 能 力 模型 


2.1 目标 函数 
航路 网 络 短期 流量 管理 的 目标 是 使 整个 航线 网 络 的 航空 
器 总 航行 费用 最 小 ， 故 可 将 目标 函数 设置 为 


min:Z=>》> xc OOD) (0) 


中 : zZ 表示 系统 总 航行 时 间 , 目标 为 求 航路 网 络 各 条 航路 费 

之 和 最 小 ， 也 即 整个 航线 网 络 总 航行 时 间 最 小 ， 其 中 随时 

与 空中 交通 流量 变化 的 cj.D) 可 由 式 (2) 求 得 。 

2 约束 条 件 
在 固定 的 安全 间隔 下 , 航路 rs 每 次 只 能 有 1 股 交通 流通 

过 ， 设 航路 网 有 nn 个 节点 ， 对 于 Yrs=12… 可 得 


(5) 


| I 


D2, n=1 (7) 
每 一 股 交通 流 在 第 i 个 时 间 段 又 只 能 通过 一 条 航路 , ， 
Vf'=1,2,..,n 可 得 ， 
bp =] (8) 


在 将 交通 流 指派 到 各 个 航路 上 时 ， 需 考虑 航路 的 容量 和 
即将 到 达 的 下 一 机 场 节点 的 容量 , 设 G8 表示 航路 rs 的 容量 ， 
G; 表示 机 场 节点 s 的 容量 ， 故 可 添加 容量 约束 ， 


之 二 六 (D< Cg » 


2 hE (10) 


关于 机 场 容量 与 航路 容量 的 确定 方法 ， 在 很 多 文献 中 已 
给 予 介 绍 和 说 明山， 本 文 不 再 资 述 。 


3 ”延误 优化 模型 的 改进 匈牙利 算法 


3.1 经 典 匈 牙 利 算 法 的 不 足 
经 典 匈牙利 算法 能 够 有 效 地 求解 一 些 工程 应 用 问题 ， 
体 步 又 可 见 文献 [15]。 但 经 典 算法 对 于 ATFM 问题 具有 以 下 
不 足 : 
a) 空 中 航路 网 是 一 个 复杂 的 巨 系 统 ， 各 航路 上 的 费用 转 
换 为 效率 矩阵 后 会 使 得 计算 的 时 间 复 杂 度 很 高 ， 从 而 不 能 应 
用 于 实际 的 管制 运行 中 。 
b) 随 着 空中 交通 流量 的 增 大 ， 航 路 网 络 会 出 现 拥 挤 ， 各 
条 航路 的 效率 值 〈 拥 挤 程 度 ) 会 出 现 变化 ， 经 典 匈牙利 法 将 
效率 矩阵 视 为 常数 对 于 空中 航路 网 的 流量 指派 是 不 合理 的 ， 
需要 引入 上 文中 提 到 的 效率 函数 。 
3.2 算法 的 改进 
1) 算 法 改进 的 思想 
针对 3.1 节 中 的 第 一 个 不 足 ， 改 进 的 思想 是 尽 可 能 减少 
算法 的 迭代 次 数 在 第 二 步 加 圈 的 零 元 素 个 数 mw 与 效率 矩阵 


的 阶 数 > 相等 的 概率 最 大 。 所 以 可 以 作出 一 些 改善 ， 如 下 文 
中 的 第 a)b) 步 。 

针对 3.1 节 中 的 第 二 个 不 足 ， 改 进 的 思想 是 用 上 文中 定 
义 的 费用 函数 描述 矩阵 ， 即 效率 矩阵 中 的 值 是 可 变 的 ， 随 着 
流量 和 时 间 的 改变 而 变化 。 


2) 改进 后 算法 的 详细 步 又 

基于 以 上 思路 ， 针 对 空中 航路 网 结构 复杂 、 交 通 流 时 变 
后 各 航路 通行 能 力 出 现 变化 的 特点 ， 改 进 的 匈牙利 算法 步 又 
如 下 所 示 。 
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a) 初始 化 。 令 航路 网 路 各 条 航路 上 的 交通 流量 流量 也 出 现 增 加 。 以 双流 机 场 和 江北 机 场 两 个 繁忙 机 场 为 例 ， 


不 =0,vr,seE ， 从 而 求 得 航路 网 路 的 初始 阻抗 2* ; 令 时 间 段 i ”算法 能 将 进 场 通行 能 力 提 升 8.372%， 离 场 通行 能 力 提升 
需 分 配 的 流量 久 ， 和 迭代 次 数 二 1 。 8.999% 。 
b) 效 率 矩 阵 的 预 处 理 。 比 较 效 率 和 矩阵 每 行 最 小 元 素 的 个 表 2 优化 前 的 离 场 交通 量 
数 Nine 与 每 列 最 小 元 素 的 个 数 Ny 8 在 Ne > Nw 3 效率 矩阵 Table 2 ”Number of departure aircrafts before optimization 
从 列 开始 迭代 ， 若 Ni <WN,s, ， 效 率 矩 阵 从 行 开 始 迭 代 。 成 都 双流 重庆 江北 桂林 两 江 南宁 吴 圩 张家界 荷花 贵阳 龙洞 堡 
0c) 流量 分 配 。 按照 经 典 的 匈牙利 算法 的 第 三 步 ~ 第 五 步 对 SO 2 9 3 0 
流量 进行 指派 ,得 到 时 间 段 ;的 流量 分 配 结果 若 在 航路 -ss 上。 7 2 , | 1 
分 配 得 流量 和 学 及 > 且 之 二 广 0>Ci ， 则 未 分 配 的 剩余 流 T0020 2 5 8 7 2 1's 
四 2:00-13:00 27 25 9 2 3 14 
量 〈 分 配 后 剩余 的 流量 ) 进入 下 一 次 迭代 ; 若 未 达到 容量 约 3:00-14:00 ”25 30 10 9 3 19 
束 ， 该 段 航路 仍 有 剩余 容量 。 4:00-15:00 20 26 9 0 3 16 
d) 更 新 航路 阻抗 。 将 分 配 得 到 的 流量 人 ,vrssE 加 载 到 航 7 2 ， 
路 网 络 , 得 到 新 的 航路 阻抗 (7.(0) , 并 转换 成 新 的 效率 矩阵 。 oso 2% 
e 收 敛 判 断 。 设 fs 为 各 个 时 间 段 需 分 配 的 总 流量 ， 若 8:00-19:00 24 25 6 二 | 16 
FE， 0 
21:00-22:00 -33 35 6 6 4 19 
4 ”仿真 算 例 22:00-23:00 ”35 30 4 2 4 17 
23:00-23:59 40 35 3 7 2 11 
4.1 算 例 构造 ee | 表 3 优化 后 的 进 场 交 通 量 
、 中 南 和 西南 空 管 本 所 条 和 最 繁忙 的 空域 Table 3 Number of approach aircrafts after optimization 
故 选取 该 地 区 空域 验证 本 文 的 模型 与 算法 。 算 例 数据 来 源 于 成 都 双流 重庆 江北 桂林 两 江 南宁 吴 二 张家界 荷花 贵阳 龙 泣 保 
中 南 和 西南 空 管 局 提供 的 2016 年 11 月 空域 流量 统计 数据 ， 人 号 
以 MATLAB2017a 作为 运算 平台 ， 验 证 上 文 的 模型 与 算法 。 G00 i000 2 25 ]2 各 i 20 
4.2 仿真 结果 10:00-11:00 19 28 9 16 1 35 
4.2.1 通行 能 力 的 增加 11:00-12:00 20 24 4 15 1 5 
根据 2016 年 的 统计 数据 , 剔除 诸如 军 航 活 动 、 和 危险 天 气 12:00-13:00 24 29 11 6 3 8 
等 流量 异常 数据 ， 仿 真空 域 六 个 主要 民用 机 场 : 成 都 双流 机 0 和 20 0 4 
场 、 重 庆 江北 机 场 、 桂 林 两 江 机 场 、 南 宁 吴 圩 机 场 、 张 家 界 0 5 1 3 3 3 
1 


荷花 机 场 和 贵阳 龙洞 堡 机 场 的 日 均 进 场 和 离 场 通行 能 力 , 见 。 “3 3 


,| 、 、 王 Mis 16:00-17:00 17 24 3 11 1 8 
表 1、2。 由 表 1、2 和 历史 数据 可 得 成 都 双流 机 场 和 重庆 ee ; , 
北 机 场面 临 较 大 的 管制 运行 压力 。 Re 
1 优 化 的 . 场 ; 有 18:00-19:00 23 24 9 于 3 3 
1 \ 入 ] 用 
个 Wo 19:00-20:00 24 23 10 13 2 3 
Table 1 Number of approach aircrafts before optimization 20.00.01:00 24 3 6 9 六 9 
成 都 双流 ”重庆 江北 ”桂林 两 江 ”南宁 吴 圩 ”张家界 荷花 ”贵阳 龙洞 保 21:00-22:00 20 12 5 3 4 7 
8:00-9:00 33 27 5 8 1 19 2200 8 该 5 5 6 8 
9:00-10:00 22 22 10 8 1 17 .00 4 a 1 1 1 3 
10:00-11:00 16 27 8 16 31 
表 4 优化 后 的 离 场 交通 量 
11:00-12:00 19 21 4 14 1 12 ee 
Table 4 Number of departure aircrafts after optimization 
12:00-13:00 21 27 10 6 3 17 
成 都 双流 ”重庆 江北 ”桂林 两 江 ”南宁 吴 圩 ”张家界 荷花 ”贵阳 龙洞 保 
13:00-14:00 20 18 6 7 4 12 
8:00- 9:00 11 20 11 5 | 20 
14:00-15:00 27 23 11 12 3 13 
9:00- 10:00 23 26 yg 19 1 18 
15:00-16:00 29 21 10 9 2 18 
10:00-11:00 27 办 6 13 1 18 
16:00-17:00 15 23 3 10 1 17 
11:00-12:00 26 35 8 9 2 18 
17:00-18:00 24 22 6 8 1 13 
12:00-13:00 27 25 9 记 3 16 
18:00-19:00 21 23 8 12 3 13 
13:00-14:00 25 30 10 11 3 19 
19:00-20:00 22 22 8 13 六 13 
14:00-15:00 20 26 9 10 3 16 
20:00-21:00 22 22 6 8 2 18 
15:00-16:00 21 28 3 9 1 12 
21:00-22:00 17 12 5 3 4 7 
16:00-17:00 22 25 8 12 2 18 
22:00-23:00 8 12 5 5 6 8 
17:00-18:00 24 26 8 15 1 16 
23:00-23:59 4 3 1 1 1 3 
- 一 ee ee 18:00-19:00 24 25 6 11 3 17 
使 用 改进 的 匈牙利 算法 对 上 述 6 个 机 场 进行 运算 求解 ，  ，。 ， ， 。 
可 算得 表 3 和 4 所 示 的 优化 后 的 通行 能 力 ， 两 个 空域 繁忙 的 jw » » ;， ， 
机 场 : 成 都 双流 与 重庆 江北 的 通行 能 力 有 明显 的 提升 。 但 是 2100.2200 33 35 6 6 4 21 
于 桂林 两 江 机 场 、 南 宁 吴 圩 、 张 家 界 和 荷花 和 贵阳 龙洞 堡 虽 22:00-23:00 35 30 4 12 4 到 
然 本 身 不 存在 空域 拥挤 ， 但 整个 空域 网 络 是 一 个 系统 ， 拥 挤 23:00-23:59 40 35 3 2 11 
空域 的 优化 使 更 多 的 航空 器 避免 了 地 面 等 待 ， 这 4 个 机 场 的 
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使 用 变化 的 费用 函数 替代 效率 矩阵 不 变 的 系数 。 


4.2.2 延误 时 间 的 优 


的 两 个 机 场 优化 前 


双流 机 场 的 延误 水 平平 均 下 降 3 min， 
延误 下 降 2 min。 国 


n 


构成 延误 ， 算 法 将 高 峰 时 间 段 每 架 飞 必 
在 13 min 以 内 。 


通过 式 (4) (5) 计算 航路 网 络 的 延误 ,仿真 空域 最 繁忙 以 西南 空 管 局 与 中 南 空 管 局 所 辖 部 分 空域 作为 算 例 ， 以 
与 优化 后 的 延误 水 平 如 图 1、2 所 示 。 成 都 。 历史 数据 计算 仿真 空域 上 主要 机 场 的 通行 能 力 ， 并 使 用 建立 

重庆 江北 机 场 的 平均 ”的 模型 与 改进 的 匈牙利 算法 对 所 辖 空域 进行 优化 ， 仿 真 结果 

际 民航 组 织 定 义 延 迟 起 飞 超 过 15 min 便 “表明 ， 优 化 模型 与 改进 的 匈牙利 算法 可 以 有 效 地 提升 航路 网 

用 的 平均 延误 大 多 控制 。 络 通行 能 力 ， 降 低 每 架 航 空 器 的 平均 延误 且 算 法 相 比 于 经 典 


匈 


牙 利 算法 具有 效率 优势 ， 更 适用 于 实际 的 管制 运行 。 
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